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Abstract 

The results presented at the 28th International Cosmic Ray Conference 
on neutrino and muon physics are summarized. Neutrinos and muons provide 
a huge amount of information on particle interactions up to very high energies 
and on fundamental particle properties. Results on neutrino oscillations in the 
atmospheric and solar u sectors are summarized. Oscillations are well established 
in both sectors, and a more precise determination of oscillation parameters is 
requested in the next future. 

Neutrino telescopes taking data and under construction presented numer- 
ous results. Neutrinos as probes of the Universe are hopefully going to open, 
together with gravitational waves, a new era for Astrophysics. Cosmology has 
entered the precision era and the Dark Matter quest is still an open problem. 
Direct and indirect searches are complementary approaches to the problem. 

The results presented at this conference confirm that Astroparticle Physics 
and, in particular. Neutrino Physics are leading fields in fundamental research. 

1. Introduction 



The High Energy Phenomena Sessions HE 2.1-5 of this conference are 
devoted to: muon and neutrino experiments and calculations, neutrino telescopes 
and new projects. The Sessions HE 3.2-4 are dedicated to more exotic searches, 
such as dark matter and new particles, proton decay and cosmology, both from 
the theoretical and the experimental point of view. The total amount of 139 
talks and posters cannot be entirely summarized here. In the attempt to provide 
an organic overview on the subjects of these sessions, some of the works are not 
mentioned. I do apologize with the authors. Whenever possible experimental 
results and calculations are compared. 

Sec.|2I]is devoted to neutrinos and muons of atmospheric origin. In Sec. 12. 1.1 
the current status on our experimental and theoretical knowledge of the atmo- 
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spheric neutrino beam and its impact on the oscillation scenario is discussed. 
Super-Kamiokande (SK) has presented results on atmospheric neutrinos (HE 2.2) 
and on z/^ — * Vr oscillations. In Sec. 12. 2. | the status of atmospheric neutrino 
calculations is summarized (HE 2.4), with particular attention to the primary 
cosmic ray (CR) flux, a crucial input to atmospheric cascade calculations. 

Muon flux measurements (HE 2.1), currently, are mainly motivated to 
benchmark atmospheric shower development codes and hence atmospheric v cal- 
culations. In fact muons and neutrinos come from the same decay chains. In 
Sec 12.3.1 the results presented at the conference on muon fluxes and charge ratios 
are summarized. 

Sec. ED is devoted to neutrinos of astrophysical origin, not produced in the 
Earth atmosphere. Among these, solar neutrinos are used to investigate neutrino 
properties as well as the Sun itself. In Sec. 13.1.1 results from SNO and SK (HE 
2.2) are summarized. The solar neutrino experiment results together with the 
KamLAND long baseline reactor experiment results are providing a convincing 
solution to the solar neutrino problem on z/g disappearance. 

Currently solar neutrinos and events registered by Kamiokande and 1MB 
[1] a few hours before the optical identification of SN1987A are the only astro- 
physical neutrinos detected so far. In Sec. 13.2.1 limits on supernova (SN) collapse 
are summarized. On the other hand, no neutrino of astrophysical origin has yet 
been detected above the GeV scale in the background of atmospheric neutrinos. 
Nevertheless hopefully in a few years we will enter the neutrino astronomy era, 
thanks to the numerous efforts on the construction of huge neutrino telescopes to 
which Sec. 13.3.1 is devoted (HE 2.3). Results from various neutrino telescopes on 
point-like and diffuse sources are presented, including sensitivities expected for 
experiments under construction and R&D. 

In Sec. 1^ some of the results and sensitivities on dark matter searches 
through the detection of gammas, anti-protons, positrons and neutrinos by satel- 
lites, ground-based arrays and neutrino telescopes will be summarized. 

2. Atmospheric Neutrinos and Muons 

2.1. Atmospheric Neutrinos and Super-Kamiokande 

Super-Kamiokande is an ultra-pure water Cherenkov detector, with a fidu- 
cial volume of 22.5 ktons and a 40% phototube (PMT) coverage (11,146 51 cm 
Hamamatsu PMTs looking to the inner detector and 1885 20 cm PMTs of the 
outer veto) [2]. A reaction chain in Nov. 2001 destroyed 50% of the PMTs. This 
unfortunate event concluded the SK-I phase with a total amount of 1489 days 
of data taking (91.7 kt yr). In Dec. 2002 the reconstruction was completed and 
since Jan. 2003 the experiment has been taking new data and the K2K neutrino 
beam has been on. 

Even though no experiment has yet measured the oscillatory pattern. 
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Table 1. Statistics and flavor ratios for the Sub-GeV and Multi-GeV samples. Mea- 
sured and expected fluxes (for no oscillations) are given for upward through-going 
/is and stopping ^s 



Sample 


e-like 


/i-like 


(^J./e)MC 


Sub-GeV 


3353 


3227 


0.649 ± O.OWstat ± 0.051,j,, 


Multi-GeV+PC 


746 


1564 


0.700 ±^omo ±0.083,^,, 


Sample 


Data 


Measured flux 


Expected Flux [5] 




(1645.9 d) 


10"^^ cm"^ sr~^ 


10~^^ cm~^ sr~^ 


Stopping fis 


463 


0.41 ± 0.02stat ± 0.02,^, 


0.61 ±0.14 


Through-going /is 


1843 


1.70 ±0.04,t„t± 0.02,^, 


1.57 ±0.35 



which would unequivocally establish an oscillation phenomenon, the strength of 
Super-Kamiokande result relies on the measurement of topologies belonging to 
different energy ranges and of both electron and muon flavors. The updated 
analysis presented at the conference introduced some reflnements in the event 
reconstruction and in particle identiflcation (on which the systematic error is 
at the level of 1% as estimated from Monte Carlo studies, from a test of the 
KEK proton synchrotron [3], and conflrmed by the near 1 kton detector of K2K 
[4]). Also improvements on the upward through-going muon selection have been 
adopted. Moreover, the Monte Carlo (MC) generator has been updated using the 
calculation by Honda et al. [5], which adopts the flt to primary CR measure- 
ments presented at ICRC2001 [6]. Also the quasi-elastic and Itt production cross 
sections were improved thanks to the K2K data. 

Tab. 1 summarizes the statistics and the values of the flavor ratio for the 
Sub-GeV fully contained (FC) sample {E^is < 1.33 GeV) and the Multi-GeV (FC 
with E^is > 1.33 GeV) and partially contained (PC) sample (< E^, >~ 10 GeV). 
For the Multi-GeV sample the up/down asymmetry for /x-like events deviates 

from zero (expected in the no oscillation case) by 9.5 a: A^_iike = j = 

—0.289 ± 0.028stat ± O.OOisys- On the other hand, for e-like events Ae-uke = 
—0.020 ± 0.043stat ± 0.005sys. These measurements provide a robust evidence 
in favor of Ur oscillations, as will be discussed in Sec. 12.2.1 Muon neutrino 

charged current (CC) interactions in the rock below the detector produce stopping 
muons (< El, >~ 10 GeV) and upward through-going muons (< E^, >~ 100 
GeV). The zenith angle distribution of these events shows a distortion compatible 
to oscillations. Tab. 1 summarizes the statistics, the measured and expected fluxes 
also for these samples. 

The SK-I updated analysis flrstly presented at ICRC2003 [2] leads to a 
preliminary best flt Am^ value of 2 ■ 10"^ eV^, at maximal mixing, about 20% 
smaller compared to the past result of 2.5 • 10"'^ eV^ [7]. Allowed regions are 
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compared in Fig. 1, where the final result from MACRO [8] and Soudan 2 [9] 
are presented too. Also the Kamiokande [7] result is shown: the lowering of SK 
region increases the discrepancy between the SK and its precursor result which 
could be clarified through a reanalysis of the Kamiokande data. The 90% c.l. 
region is sm 29 > 0.9 and 1.3 < Am^ < 3.0 • 10"^ eV^ while previously it 
was sm'^2e > 0.9 and 1.6 < Am^ < 3.9 ■ 10"^ eV^ The value of Am^ has 
an impact on the expectations of Long Baseline Experiments, particularly on 
the CNGS appearance experiments, OPERA and ICARUS. Considering a beam 
intensity larger than a factor 1.5 compared to the nominal beam, the CC Ur events 
expected in ICARUS T3000 (2.35 kton active mass) in 5 yrs for the decay modes 
T ^ e , p vary from ~ 7 to 18 between Am^ ~ 1.6 — 2.5 ■ 10"'^ eV^ and maximal 
mixing [10], while the background is ~ 1 event. 

Since one of the changes is related to the p fiux adopted by the experiment, 
the problem, addressed in Sec. 12. 2. [ on how well we do know the atmospheric u 
fiux in the various energy ranges arises, even though it is clear that also the 
analysis refinements must have a relevant role. 

In the low energy range (PC events) SK investigated the geomagnetic field 
effects on u fiuxes through the measurement of the azimuth distribution affected 
by the East- West anisotropy in the primary CR fiux. The anisotropy results in 
larger (smaller) fiuxes of z/'s traveling toward East (West) [11]. This offers the 
possibility to validate calculations at < 2 GeV. SK is also sensitive to the finer 
effect of an enhancement of the asymmetry for e-like events with respect to //-like 
events that can be correctly taken into account by 3-D MC calculations {u's are 
not collinear to their parents) [12]. 

SK can provide information on the oscillation channel: z/^ T-'sterUe oscil- 
lations are disfavored at 99 % c.l. through the non observation of neutral current 
(NC) suppression using the vr" sample which now is measured with a reduced 
systematic error (from 30% to 9%) thanks to the K2K measurement of = ^ 
[13]. Moreover some indication on r appearance can be extracted on a statistical 
basis [14]. 

2.2. The Knowledge of the Atmospheric Neutrino Beam 

Exploiting different experimental techniques, SK, MACRO and Soudan 2 
provide a robust evidence in favor of t'^ Ur oscillations, since the flavor ratio, 
the up-down asymmetry and the zenith angle distributions are robust quantities 
against theoretical errors on the atmospheric u flux calculation, known at the 
level of < 5%. As a matter of fact, at energies below about 2 GeV, when all 
muons decay in the atmosphere, 2 z/^ and 1 are produced in vr and fi decay 
chains. Moreover, far from regions where geomagnetic fleld effects are relevant, 
that is for E,y >2 GeV, the u flux is up-down symmetric {^iy{E, 9) = ^u{E, tt — 
6), where 6 is the zenith angle), due to the atmosphere and Earth spherical 
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Fig. 1. The 90% c.l. updated (dashed red Une) and the previously (dotted thin red 
hne) allowed regions from the combined analysis of all samples for SK-I data. The 
Kamiokande (black solid thin line) , MACRO (black solid line) [8] , Soudan 2 [9] (blue 
'irregular' solid line) are shown and indicated by names and arrows. SK (MACRO) 
updated best fit is indicated by a red circle (black star) at Am^ = 2 • 10~^ eV^ 
(2.3 • 10"'^ eV^) and maximal mixing. 



symmetry. For higher energy events, i.e. the through-going muons, the shape 
of the angular distribution obtained using different calculations agrees at the 
level of ~ 5%, where the remaining uncertainty is due to the knowledge of the 
competition between decay and interaction for pions and kaons, seasonal effects 
and variations in atmospheric profiles, and to the uncertainty in the spectral 
slope of the primary flux. On the other hand, using the comparison between 
different atmospheric neutrino calculations, it can be estimated that the absolute 
normalization of fluxes is still affected by 15%, 30% errors, respectively in the < 10 
GeV and > 100 GeV energy ranges. This uncertainty is large even if calculations 
are benchmarked against atmospheric muons at sea level and in the atmosphere 
(see Sec 12. 3. j) . The main errors come from the primary cosmic ray spectrum 
and the hadronic interaction models. During the last years, authors of different 
MC computations started a comparison work between interaction models which 
led to code updates. For instance, the Bartol group presented at the conference 
the updated generator TARGET [16, 17]. Other groups use more sophisticated 
interaction and transport codes, such as FLUKA [18], that has been extensively 
benchmarked on accelerator data and muons (see for instance the comparison of 
meson multiplicity distributions in xiab for generators used in CORSIKA [19] with 
p-^Be, p-^'^N data), the Japanese group updated the interaction model adopting 
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DPMJET3 [5]. Wentz et al. adopted some of the CORSIKA generators, among 
which DPMJET II. 5 [20], while Liu et al. adapted a parametrized hadronic model 
pubhshed in 1989 [21]. It should be understood why Liu et al. calculation seems 
to predict a larger East- West asymmetry for /i-like events than e-like, contrarily 
to expectations (see Sec. I2.1.|) . All of these codes are now 3-dimensional: the 
introduction of transverse momenta produces a considerable enhancement at the 
horizon for energies < 2 GeV, first found in [18]. A review on atmospheric neutrino 
calculations is in [15]. 

The other important input in atmospheric neutrino calculations is the 
primary CR spectrum. The uncertainties affecting the CR spectrum are energy 
dependent. The recent fit presented at ICRC2001 [6] was an attempt to unify 
this input between different calculations. At energies < 200 GeV/nucleon it relies 
on the recent measurements of the AMS-01 Space Shuttle flight [22] and of the 
balloon experiment BESS 98 [23]: for protons they are in agreement within 5%, 
while for He there is some discrepancy at the level of 15%, which at these energies 
reflects on a 3% uncertainty on the all-nucleon flux [15]. Nevertheless, it remains 
to be understood the disagreement of the proton measurements by the CAPRICE 
balloon flights (1994 [24] and 1998 [25]) which are about 20% lower than AMS-01 
and BESS 98 results. 

At high energy (> 10^ GeV/nucleon) the JACEE [26] and RUNJOB [27,28] 
data have large errors leaving room for large variations of the fltted spectral slope. 
While at energies < 100 GeV/nucleon, the ICRC2001 flt [6] changed by less than 
10% compared to the previous Bartol 96 flt used in [29], at higher energies the 
ICRC2001 flt [6] is steeper (~ E"^-^^) than the Bartol 96 flt (~ E"^-'^^). In 
Fig. 2 (on the left) the data and both flts are shown. ATIC has presented at this 
conference preliminary data [30] from the long duration balloon flight that flU the 
gap in energy between previous experiments. These data are crucial to understand 
the slope of the spectrum and also if protons and He have the same slope, as it 
is discussed in [31]. Preliminary results indicate a i?"^'^^ preference and that He 
spectrum has almost the same slope of the proton one. A preference for harder 
spectra for He (closer to JACEE result than to RUNJOB) is also found using 
EAS-TOP/MACRO and the analysis of coincidences between the 2 experiments 
[32]. 

The effect on u fluxes of changing the CR flux slope above 100 GeV/nucleon 
from E~'^''^^ to E~'^''^^ has been investigated by the Japanese [5] and FLUKA [18] 
groups. Using both flts with the same interaction code, in the case of FLUKA, 
the ICRC2001 flt produces + fluxes (averaged over the entire hemisphere) 
never larger by 10% between 0.1-1 GeV, lower by 5% at 10 GeV and by 20-30% 
between 100 GeV-1 TeV than the flux with the Bartol 96 flt. 

One more indication on the fact that the ICRC2001 flt could be too steep 
comes from SK and MACRO through-going /i data. As a matter of fact, if the 
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measured SK zenith angle distribution of through-going muons is fitted letting 
the Honda et al. flux [5] normalization free, the best fit is obtained for a normal- 
ization value larger by ~ 25% than what predicted by [5] using the ICRC2001 fit. 
A similar result is obtained using MACRO upward through-going muons [8,18]. 
Nevertheless it is hard to understand the effect of the energy dependent uncer- 
tainty on the normalization in the evaluation of the best fit parameters. In fact 
the Super-Kamiokande fitting procedure contains various parameters related to 
the normalization and the slope of the v spectrum [2]. In Fig. 2 (on the right) 
MACRO through-going muons (821 events selected for oscillation studies) [8,18] 
are compared to the FLUKA fluxes with the ICRC2001 and 96 fits, the Bartol 
96 flux and the Honda et al. flux [5] using the ICRC2001 fit. Clearly the data 
seem to prefer the fluxes using the 96 fit. From the plot it is also noticeable the 
good agreement between the calculations using two different interaction models 
(FLUKA and DPMJET3) and the same CR spectrum (ICRC2001 fit). 

2.3. Muon Data at Sea Level and in the Atmosphere 

Many balloon and ground based experiments presented measurements of 
muon fluxes at sea level and in the atmosphere. A summary of sea level data, 
including previous measurements is presented in Fig. 3 in terms of vertical muon 
flux and charge ratio /i^/yU", compared to the world average band estimated 
in [33]. At energies > 5 GeV muon data are not affected by solar modulation 
and geomagnetic cut-off effects, while for lower energies these effects must be 
considered, as shown for instance by the comparison between the data taken at 
Ft. Sumner by CAPRICE 97 [34] and BESS-01 [35]. 

Preliminary results on atmospheric muons from LEP experiments have 
been presented at the conference. The L3-|-Cosmics [36], which uses the L3 mag- 
netic fi spectrometer and a scintillator setup to provide the fi arrival time, has 
measured the muon spectrum in a wide energy range (20-2000 GeV) with a sys- 
tematic error of 2.6% at ~ 100 GeV (mainly due to the uncertainty on the knowl- 
edge of the 30 m thick molasse overburden) increasing to 15% at 2 TeV due to 
the momentum resolution. Using ALEPH hadronic calorimeter and TPC, the 
Cosmo- ALEPH experiment [37] presented the muon flux and charged ratio mea- 
surement in the range 70-2500 GeV. The results are preliminary since the flux is 
normalized to the world average and the experiment is investigating trigger effi- 
ciencies. Both results are of great interest for benchmarking atmospheric shower 
development codes (see Sec. 12. 2. |) up to high energies. The L3+C measurement 
has already been compared to the Bartol calculation in [16]. Other comparisons 
between calculations and balloon flight ascent data at < 20 GeV have been shown 
in [41,42,21]. The interest of data taken at the top atmosphere (see for instance 
[41]) relies on the possibility of testing the first interactions, when the shower has 
not yet developed. 
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Fig. 2. On the left: The 3 upper black lines are the all nucleon primary spectrum 
resulting from the ICRC2001 fit [6]. The middle blue group of 3 lines are the proton 
primary flux and the lower group the He flux. The upper and lower lines of each 
group are the maximum and minimum fluxes from the errors of the fit, respectively. 
The dashed line is the all nucleon flux adopted by the Bartol group in 1996 [15]. 
The data points for protons and He are from JACEE [26] and RUNJOB (open 
symbols: [27], full symbols: [28]), AMS-01 [22], BESS 98 [23], CAPRICE 98 [25], 
CAPRICE 94 [24]. On the right: upward through-going muon flux {E^ > 1 
GeV) measured by MACRO [8,18] (red symbols) compared to the Bartol 96 [29] 
(solid green line), to the FLUKA calculation using the Bartol 96 flux (black dotted 
line) and the ICRC2001 fit (pink dashed line), to the Honda et al. flux using the 
ICRC2001 fit [5] (blue dash-dotted line). Predicted curves include oscillations with 
maximal mixing and Am2 = 2.3 • lO^^ eV^. 

3. Neutrino of Astrophysical Origin 

In this section results concerning measurements and searches for neutrinos 
not produced in the Earth atmosphere, such as solar neutrinos and neutrinos from 
cosmic sources, are summarized. The only exception is for the reactor experiment 
KamLAND, that is included in the solar u session due to the impact of its results 
on the solar u problem. 

3.1. Solar Neutrinos 

In writing this section on solar neutrinos it cannot be ignored the fact 
that after the conference the SNO collaboration has published results on the salt 
phase [43] and numerous interpretations have already appeared. SNO, a 12 m 
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Fig. 3. On the left: Vertical atmospheric muon flux at sea level (altitude corrected 
using the formula in [33]) vs momentum. The yellow band is the flux world average 
in [33]. Data: L3+Cosmics (black triangles) [36], red full circles (Cosmo- ALEPH 
[37]), Rastin 84 (ref. in [33], green triangles with the highest normalization at 100 
GeV). At energies < 10^ GeV: BESS-01 at Ft. Sumner [35] (smah red triangles), 
BESS 95 at Tsukuba (small blue reverse triangles) and BESS 97-99 at Lynn Lake 
(green crosses) [38], BESS-99 at mountain altitude at Mt. Norikura (cyan squares) 
[39], CAPRICE 94 (pink reverse triangles not distinguishable from BESS 95) and 
CAPRICE 97 (black small circles) [34]. On the right: Muon charge ratio: same 
symbols as the plot on the left, including also OKAYAMA (large pink squares) [40]. 



diameter acrylic vessel filled with 1 kton pure D2O seen by 9500 PMTs, exploits 
the solar v detection above a ~ 5 MeV threshold through 3 reactions: neutrino 
elastic scattering (ES), mainly sensitive to z/g since a^^ ~ ^^vr,v^ CC and all 
flavor V NC interactions on deuterium. The ratio of CC to NC interactions 
provides the electron survival probability. The results of the first phase of the 
experiment, also combined with the SK precise measurement of ES showed that 
about 2/3 of i^e convert into another active flavor or i^^, indicating hence a clear 
preference for MSW effects [44] and a 5.3(T rejection of the null hypothesis on flavor 
transformation. After increasing sensitivity to NC with the addition of ~ 2 tons 
of salt, the experiment is now able to claim evidence for MSW effects in solar 
neutrinos. These results represent a fundamental step forward in understanding 
solar V oscillations through the precise determination of the total active (and 
hep) flux in agreement with [45]. 

SK [46] presented the updated analysis of 1496 days data taking with 
Eth = 5 MeV. In this period 22385 solar ps were detected through the recoiling 
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electron in ES interactions with sensitivity to v energy and direction. The day- 
night asymmetry parameter is (J^^^-^j^ = — 0.021 ±0.020±g;o}25 which is consistent 
with zero, and the electron energy spectrum is consistent with being fiat with 44% 
C.I.. No significant time variation or energy distortion appears in the data favoring 
LMA (large mixing angle) solution. Fig. 4 (on the left) shows the allowed region 
at 95% c.l. including SK information on rate, spectrum and time-variations, 
SNO measurements of CC and NC available at the time of the conference [44], 
radiochemical experiments and the region allowed by the measurements of the rate 
and of the spectrum of reactor z/s in KamLAND [47]. About 80% of the muon 
neutrino fiux reaches the KamLAND detector (1 kton scintillator with a 34% 
PMT coverage) from reactors at 140-210 km of distance with < >~ 3 MeV. 
This disappearance experiment observed a deficit of events with respect to the no- 
oscillation expectation of = 0.611 ± O.OSS^fat ± 0.041^^^ and measured 
the energy spectrum of prompt positrons from the reaction v^. ^ "p ^ ^ n. 
The LMA region is divided into two allowed parts including KamLAND results 
(if CPT is conserved). After the results of SNO salt phase [43] the higher Am^ 
region disappears at 99% c.l. and the best fit point is Am^ = 7.1 • 10~^ eV^ 
and tar?Q = 0.41. KamLAND [47] unfortunately has not presented new results 
since most of the Japanese reactors have been running at much reduced efficiency 
due to safety controls. Preliminary results on anti-neutrinos produced in the 
Earth crust in U/Th decays (geo-neutrinos with between 0.9-2.5 MeV) were 
presented posing a limit on the Earth heat source of < 110 TW (95% c.l.), strongly 
dependent on the element distribution particularly of the Japanese Island Arc. 

A search for possible time modulations has been presented by SK after 
some papers appeared claiming for periodicities very close to the time interval in 
which published data are binned (10 days) [48]. No significant periodicity (except 
for the long term modulation due to the eccentricity of the Earth orbit around 
the Sun) is found and modulations of 10-100 (30-100) days with amplitude larger 
than 10% are excluded at 95% c.l. binning data in 5 (10) days. 

Even though the LMA solution seems well established, the current pre- 
cision on parameters still allows for sub-dominant processes. Limits on solar z/g 
appearance due to the combined effect of spin flavor precession in presence of a 
large enough solar magnetic fleld and neutrino magnetic moment * (evolving z/g 
into z/^ , Ur) and of flavor oscillations (converting ,Ur ^e) were presented by 
SNO [50] and SK [51]. This search addresses the fundamental question on the 
nature of neutrinos, Majorana or Dirac particles. As a matter of fact, z>eS would 
be detected only if z/s are Majorana particles, since if they are Dirac ones the elec- 
tron left-handed neutrinos would convert into sterile (hence not detectable) right- 
handed ones. The reaction exploited by SK is inverse (5 decay, where positrons 

*Current best limit on the v magnetic moment is from the MUNU experiment < 10~^''/iB 
(90% c.l.) [49]. 
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Table 2. Upper limits (90% c.l.) on the flux of Ue in percentages of the Standard 
Solar Model Vf, flux [45]. Live days and positron energy ranges are indicated. 



SNO [50] 


SK [51] 


KamLAND [49] 


1.02 


0.8 


0.028 


{Eth = 5 MeV, 306 d) 


{E^in = 8 MeV, 1496 d) 


{E e [8.3, 14.8] MeV, 185.5 d) 



cannot be separated by electrons and the coincident 2.2 MeV photons from neu- 
tron capture cannot be detected since their energy is below threshold (contrarily 
to what happens in KamLAND [49]). A 93% contamination is expected from 
spallation backgrounds [51]. On the other hand cleaner signatures are found in 
SNO [50] from the reaction Uf. + d~^e~^ + n + n thanks to 3-fold and 2 fold (e"''n, 
nn) coincidences. Nevertheless, the efficiencies are only ~ 1% and 10% respec- 
tively. In Tab. 2 results are summarized including the post- conference KamLAND 
result [49]. Furthermore, it should be considered that alternative Ue sources could 
be WIMP annihilation in the Sun, relic SN neutrinos ([52], see Sec. 13. 2. |) and 
neutrino decay. 

3.2. Neutrinos from Stellar Collapse 

About 99% of the binding energy (~ 3 • 10^^ erg) of a collapsing star 
goes into neutrinos: Ue are produced during the neutronization phase which lasts 
about 10 ms, and neutrinos of all flavors during the thermalization phase of ~ 10 
s. Equipartition of energy luminosity between different u flavors and Fermi-Dirac 
spectra are expected, and average energies of < Ey^ j^^ „^ >~ 13, 16, 23 MeV, even 
though recently detailed simulations show that < Ey^ > is much closer than before 
to < Ey^^ > [53]. The largest rate of events is foreseen for inverse (3 reactions in 
water/ice and scintillator detectors. In Tab. 3 the limits on the SN rates presented 
at this conference and others are given. As a reference, the expected SN rate in 
our Galaxy is 2-4 per century. 

LVD [54] has presented an analysis on the effect on detected event rates 
due to u oscillations for the normal > 0) inverted hierarchy cases 

< 0) which could lead in conservative models to an enhancement up to 
~ 450 events for a SN at 10 kpc emitting 2.5 ■ 10^'^ erg in z/'s, or a suppression 
of events in pessimistic cases in which < E^^ > is very close to < E^^^^^ >. 
Essentially the enhancement is due to the fact that if u^ ^- oscillate into u^, since 
they have larger average energies, their spectra are harder than for Ue in the case 
of no oscillations. 
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Table 3. Limits (90% c.l.) on the SN rate in units SN/yr for various experiments. 
Distances to which hmits apply and data taking times are indicated. 



LVD [54] 


SK [55] 


AMANDA II [56] 


MACRO [57] 


0.24 


0.49 


4.3 


0.27 


< 20 kpc, 10 yrs 


< 100 kpc, 4.7 yrs 


in Galaxy, 80 hrs 


in Galaxy, 10 yrs 



3.3. The Neutrino Telescope Era 

The Neutrino Telescopes session (HE 2.3) contained numerous contribu- 
tions indicating an intense experimental activity, aiming at the detection of the 
first astrophysical i^'s with high energies. These neutrinos would constitute new 
messengers from the universe, since compared to photons which currently provide 
our best knowledge, they are less absorbed due to their weak interactions. Pioneer 
works on neutrino astronomy (see for instance the review in [58]) were initiated by 
tracking calorimeters (Kolar Gold Field) and water Cherenkov detectors (1MB, 
Kamiokande), followed by MACRO [59] and SK [60], with areas of the order of 
10^ — 10^ m^ and Eth ~ 1 GeV looking for upward through-going muons produced 
by Vfj_ interactions in the surrounding rock. 

This detection technique profits of the increase with energy of the yU. range 
and the v CC interaction. Typical astrophysical v sources are beam dumps where 
protons are accelerated on matter or gas of photons producing charged mesons, 
or decays of very massive particles, such as topological defects. At high ener- 
gies it is expected that the signal of i^'s produced in beam dumps with typical 
power law spectra from l'^* order Fermi acceleration mechanisms E'"^ — E~'^'^, 
overcomes the steeper spectrum of the atmospheric u background (~ E~^'^ above 
100 GeV). The rejection of the atmospheric /i background is achieved by measur- 
ing the up-going fi direction and by locating detectors below kilometers of matter. 
Pointing capabilities and energy measurement are relevant to reject both sources 
of background. 

The low expected event rates from astrophysical sources and current exper- 
imental upper bounds, urged the construction of huge neutrino telescopes (areas 
of the order of 0.01-1 km^, located in the South Pole ice/lake/sea water depths). 
Essentially these are 3-dimensional arrays of PMTs that allow to reconstruct /i 
tracks in water using the times of PMTs hit by the emitted Cherenkov photons 
and to estimate the energy from charge measurements. Besides cascades from 
NC, CC z/g and Ut interactions (about 1/2 of z/^ from cosmological sources os- 
cillate into Ur assuming atmospheric u oscillation parameters) are detected as 
point-like sources of light. The extreme environmental conditions of the loca- 
tions of the experiments represent a challenge and the field is rich of successes 
and drawbacks. At this conference numerous results have been presented by 
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Fig. 4. On the left: Allowed region (95% c.l.) from the combination of SK [46], 
SNO [44], Ga and CI experiments. Superimposed the darker allowed region from 
KamLAND [47]. On the right: 90% c.l. upper limits on the muon flux induced 
by neutrinos with spectrum of the form E~'^ as a function of source declination 
for SK (green circles) [60] , MACRO (red squares) [59], AMANDA-BIO (3 blue solid 
lines: the upper and lower line give the maximum variations of limits for different 
right ascension bins in the same declination band) [69], AMANDA-II (dashed line) 
foreseen limit for 1997-2000 data sample [68], ANTARES sensitivity (dotted black 
line) after 1 yr [76]. It has not been possible to apply a correction due to different 
muon average energy thresholds in this plot (SK ~ 3 GeV, MACRO ~ 1.5 GeV, 
AMANDA ~ 50 GeV). Nevertheless, the maximum of the response curves for all 
of these detectors for an E~'^ flux is -E^ ~ 1 — 10 TeV, hence the number of events 
contributing between 1-50 GeV should not make a large correction to these limits. 
Moreover, /i flux limits are less dependent on different v flux models than u limits. 

AMANDA [61] and Lake Baikal [62]. Supported by the experience acquired with 
AMANDA, the challenge at the South Pole will continue with an improvement of 
more than 2 orders of magnitude in sensitivity with the IceCube detector, whose 
construction will start in the austral summer 2004-5 and will last 6 yrs [63]. The 
community is now looking forward to see first results from experiments in the 
Mediterranean: ANTARES [64] will be completed in 2006 and common efforts 
between ANTARES, NESTOR [65] and NEMO-RD [66] could lead to the con- 
struction of a km'^ detector which would complement IceCube sky coverage. The 
need for huge detectors is also motivating R&D on new techniques alternative to 
the Cherenkov one, such is the case of the RICE radio detector [67]. 
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3.4- Status of the Experiments and Results 

AMANDA has presented first results of tlie current AMANDA-II configu- 
ration made of 677 optical modules (OM), transparent pressure resistant spheres 
containing PMTs, distributed over 19 strings inside a cylinder of 120 m diam- 
eter located at depths between 1500 and 2000 m in the polar ice. This layout 
increased the sensitivity in the high energy range of about a factor of 4 compared 
to AMANDA-BIO configuration active during 1995-7, with 302 OMs on 10 strings. 
The effective area for i/-induced muons is now almost uniform with declination 
being on average ~ 0.03 km^ at 10 TeV. The angular resolution benefits from 
the additional strings, being now about 2.3° [68] compared to about 4° [69] for 
AMANDA-BIO. The absolute pointing precision is at the level of 0.5° verified 
using the SPASE extensive air shower array (EAS) at the surface. 

Two test-beams can be used by these experiments to understand systemat- 
ics on efficiencies and acceptance: atmospheric muons and neutrinos. The shape 
of atmospheric fi vertical intensity as a function of the cosine of the zenith of 10 
hrs data is in agreement with the simulation using CORSIKA with the QGSJET 
interaction model, but the data need to be reduced by 30% to agree on the nor- 
malization [70] . This is compatible with the systematics of the experiment mainly 
due to depth- dependent ice optical properties and bubbles formed around OMs 
after drilling and hence to OM sensitivities. Also the energy spectra of both 
beams have been reconstructed with a muon energy resolution (defined as the 
standard deviation of the IoqiqE distribution) of 0.4-0.6 between 500 GeV-5 PeV. 
In particular in the 1 TeV region, AMANDA atmospheric u spectrum, having a 
fitted slope of —3.56 ± 0.20stat, agrees with the Frejus one [71]. Also a method 
has been developed to extrapolate the CR spectrum in the 1.5-200 TeV/nucleon 
region, independently of ice optical properties, OM sensitivities and detector con- 
figuration changes. This method results in a spectral index and a normalization 
for protons of -2.80 ± 0.02 and 0.106 ± 0.007 m'^ s"^ sr"^ TeV-^ respectively 
[72]. The result is in reasonable agreement with direct measurements. 

AMANDA experience is precious for the construction of the km^ scale de- 
tector, IceCube [63], with 4800 DOM (Digital Optical Modules) on 80 strings 
spaced by ~ 125 m, implemented between 1400-2400 m below the surface. DOMs 
will exploit digital readout of 10 inch PMTs storing the full waveform with a 
dynamic range of 200 photoelectrons per 15 ns. The detector will be comple- 
mented by the EAS array IceTop, 80 stations close to IceCube holes made of 2 
tanks filled with ice seen by DOMs [73]. IceCube declared effective area after 
selection cuts reaches 1 km^ at 10 TeV and the angular resolution is about 0.6° 
above 100 TeV, improving with energy. For energies > 1 PeV upgoing Ufj_ are 
non negligibly absorbed in the Earth depending on the zenith angle. Thanks to 
energy cuts getting rid of atmospheric us and /is, the experiment will be able to 
measure astrophysical down-going z/ events. 
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Results from the radio-frequency detector RICE [67] were quite impressive. 
The detector is made of 16 radio receivers with a frequency bandwidth ~ 200—500 
MHz located in holes drilled for AMANDA at depths between 100-300 m over a 
volume of about (200 m)^. The attenuation length in ice for radiation is > 1 
km. The technique exploits the detection of few ns radio pulses produced by 
V induced electro-magnetic cascades. A negative charge excess develops in the 
shower due to annihilation and extraction of e~ from the media resulting in 
a coherent Cherenkov emission (Askarian effect) proportional to the square of 
the primary particle energy. Limits comparable to AMANDA-II sensitivity were 
obtained using such a cost effective technique. Given that a better rejection of the 
anthropogenic noise could be achieved deploying receivers deeper in the ice, the 
opportunity to implement receivers also in IceCube holes should be considered. 

In the Northern hemisphere, the NT-200 Baikal detector is running at 
1100 m depth in the Siberian Lake Baikal [62], with about 192 OMs on 8 strings. 
The angular resolution is about 3°. In Mar. 2003 a prototype string has been 
deployed at about 100 m from the center of the detector with 6 couples of OMs; 
with 2 more strings the sensitivity to cascade events will be improved by a factor 
of 4. Baikal has measured 84 neutrino events (268 d for NT-200 and 70 d for 
the previous configuration NT-96) compatible with the hypothesis of being of 
atmospheric origin. 

ANTARES [64], a European project started in 1996, is deploying a 12 
string detector with a total of 900 OMs, containing 10 inch PMTs, in the Mediter- 
ranean close to the South France coasts. The effective area, estimated after se- 
lection requirements on the track reconstruction error, is > 0.02 km^ for Ey > 10 
TeV and the angular resolution achieves a limiting value of ~ 0.2° for E,y > 10 
TeV, mainly due to the transit time spread of PMTs and the diffusion of light 
in water. The 'junction box', transmitting power and data to strings, is already 
lying on the sea bed at 2500 m below the surface since Dec. 2002. It is connected 
to a 40 km long electro-optical cable deployed in Oct. 2001. Two strings have 
been deployed by a manned submarine and operated between Dec. 2002 and Jul. 
2003: one prototype string including 15 OMs with the final front-end electronics 
and the other string for environmental parameter measurements [74]. The proto- 
typing experience allowed to verify the detector design and functionality, to find 
solutions to a few occurred problems, and yielded a vast amount of measurements 
on ^^K and on a strongly variable bioluminescence activity. It was found that 
this environmental background rate is below 200 KHz for > 90% of the time. 

In Mar. 2003 NESTOR [65] deployed 10 km off of Pylos coast (Greece) 
a 12 PMT test floor of reduced dimensions compared to the hexagonal detector 
design (6 couples of up-down looking PMTs are located at 6 m from the center) 
at ~ 4000 m depth connected to shore by a 28 km cable. Cable connection 
operations performed in air on a boat and data readout and transmission were 
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successful. 

Fig. 4 (on the right) summarizes current upper hmits on //-induced /i fluxes 
for point-like sources and expected sensitivities as a function of their declination. 
The AMANDA-II sample [68] for point-like source searches consists in 699 upgoing 
events. Upper limits were given for a few sources and are at the level of a few 
10~^^ cm~^ s~^. The estimated level of limits for the 97-02 data sample is shown 
together with published results from AMANDA-BIO [69]. Also shown are the 
limits by SK [60] using 2369 through-going and stopping muons with an angular 
resolution of about 4°. The 354 showering events, used for Weakly Interacting 
Massive Particles (WIMPs) in [75], with estimated average energy of ~ 1 TeV 
could represent a 'gold' sample for astrophysics searches. The expected sensitivity 
after 1 yr of data taking of ANTARES is given in [76], where a comparison of 
binned and unbinned likelihood ratio methods is performed. 

In Fig. 5 (on the left) results on upper limits on diffuse fluxes compared 
to models are shown. The atmospheric backgrounds are rejected through energy 
cuts and n track reconstruction is required except for the Ultra-High Energy 
(UHE) AMANDA analysis [77]. The absorption of z/^ in the Earth is included. 

Cascades produced by u^., Vt CC and NC interactions are detected with 
worse angular resolution but with better energy reconstruction than events (e.g. 
for AMANDA 30° - 40° and log E resolution of 0.1-0.2 between ~ 50 TeV -100 
PeV). Cascade events are detected from the entire solid angle, the background 
discrimination being possible thanks to vertex identification. Limits on diffuse 
fluxes of all flavor neutrinos are shown in Fig. 5 (on the right). 

Many estimates on astrophysical fluxes produced after oscillations, z/^- 
regeneration processes and r propagation in the Earth have been presented up 
to the EeV energy range [81,82,83,84,85]. Even though z/t-s are not absorbed in 
the Earth due to their regeneration chain [vj. CC interaction followed by fast r 
decay which produces back and other flavor z/'s), they loose energy. Hence 
event rates are low since r leptons can be recognized at energies > 1 PeV, when 
r range in water/ice is > 50m. The r events can be unequivocally identified if 
the 2 cascades from CC interaction and decay (double bang events) are detected 
but only a few events/yr/km^ are expected [85]. Most of the UHE events (> 0.1 
EeV) can be detected from the upper hemisphere in an IceCube-like detector. For 
a GZK V model, in which UHECR make photopion interactions on the cosmic 
microwave background, ~ 50 yU and r events in a km^ per year are expected from 
the upper hemisphere and a few from the lower [81]. 

The conversion efficiency of into r for ~ 10 km of rock is much larger 
above 1 PeV than in ~ 3000 g/cm^ of horizontal atmosphere [84]. Hence event 
rates in fluorescence and Cherenkov arrays due to Earth skimming z/s (z/^-s inter- 
acting in Earth chords of the order of a few tens of kms) or mountain events (z/t-s 
crossing a few tens of kms of mountain rocks and producing a shower detectable 
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from another peak ~ 10 km far from the mountain) are more encouraging than 
event rates produced by horizontal us interacting deep in the atmosphere [84,86]. 

4. Dark Matter Searches 

After the pubhcation of WMAP data [87] cosmology entered a new pre- 
cision era. About 26.7% of the universe density is due to matter, most of which 
is dark and only 4.4% is of baryonic nature; the remaining 73.3% is due to dark 
energy. Most of the searches for dark matter (DM) presented at the conference 
are indirect in the sense that secondaries produced by dark matter annihilation 
are looked for, except for the Bulby mine experimental program on DM direct 
detection (DM particles interact in the detector) described in [88]. 

The most interesting candidates for Cold DM are WIMPs, since their 
annihilation cross section at the weak scale would account for the DM in the 
Universe. Between CDM particles, the Supersymmetric (Susy) neutralino, a lin- 
ear combination of Susy partners of the photon, Z° and neutral Higgs bosons, is 
one of the favorite candidates. Indirect searches for dark matter are performed 
by satellite, balloon and ground-based Cherenkov detectors looking for p, anti- 
deuterons, e+ and 7 excesses respect to the expected secondary fluxes produced 
by CR interactions during diffusion in the Galaxy. Two intriguing indications 
of a DM component in secondary fluxes have stimulated a lot of discussions in 
the field. EGRET has measured a diffuse gamma-ray flux in excess compared 
to standard models of primary CR interactions with the interstellar medium in- 
dicating the possible presence of a diffuse source in the Galactic Center region 
[89] . Moreover the HEAT balloon flights in 1994-95 measured a not highly signif- 
icant positron excess above secondary production models around ~ 10 GeV [90]. 
The AMS-02 experiment, to be installed on the International Space Station in 
2005, will measure all these channels up to a few hundreds of GeV [91]. As an 
example the performances for the channel are shown in Fig. 6 (on the left). 
Anti-proton fluxes from neutralino annihilation using spherical isothermal distri- 
butions of DM in the halo indicate that it would be difficult to single out the DM 
contribution in the secondary background and to constrain the Susy parameter 
space [92]. Moreover the estimates are affected by uncertainties on propagation 
models. More optimistic predictions could be obtained using different density 
profiles or hypothesis about clumpy halos. 

A cleaner DM signature with respect to excesses in diffuse fluxes can be 
monoenergetic 7 lines produced by neutralino annihilation {xx ~^ 11 ^XX ~^ 1^)- 
Searches for this signal have been presented by HEGRA [93] and Milagro [94]. 
The HEGRA system of imaging atmospheric Cherenkov telescopes surveyed the 
Andromeda galaxy (M31) looking for 500 GeV-10 TeV gamma-ray emissions with 
an energy resolution of 10%, and derived upper limits still far from models not in- 
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eluding dark matter elumpiness. Milagro water Cherenkov detector has looked for 
TeV gamma-rays from the Sun neighborhood possibly due to neutralino trapped 
in the solar system. 

Results on indirect DM searches presented by neutrino telescopes [95, 75, 
62, 102] are summarized in Fig. 6 for v induced [i fluxes produced by neutralino 
annihilation in the core of the Earth (on the right). 

5. Conclusions 

At this conference, many refinements on experimental analysis have been 
presented in the neutrino sector, even though no striking news. Atmospheric v 
experiments strongly indicate that muon neutrinos oscillate maximally into tau 
neutrinos. However, Super-Kamiokande oscillation parameter best fit value has 
changed a little after the analysis update. Naturally the question on how much 
the oscillation parameter estimate is affected by the knowledge of atmospheric v 
beam arises, a subject widely discussed at the conference where many improved 
atmospheric v calculations were presented. Still uncertainties are at the level of 
15-30% on absolute fluxes, not affecting, however, the robustness of the result in 
favor of oscillations. Solar z/s are producing striking results in these years and 
provide a useful mean to investigate fundamental properties on i/s, such as the 
question if it is a Dirac or Majorana particle. This is one of the still open problems 
which should be addressed with efforts comparable to those which are going to be 
devoted on v mixing matrix elements determination. The Dark Matter problem 
also is such a fundamental question that urge big efforts to find new evidences 
and to confirm already existing indications, such as the DAMA result [98]. 
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Fig. 5. On the left: 90% c.l. upper limits on the diffuse flux vs Ei^. Horizontal 
lines in this units are limits for E~'^ v fluxes. Experimental limits are indicated 
by the name of the experiment: AMANDA B-10 (130 days) limit up to 1 PeV [78] 
and the limit for UHE horizontal showering /is up to a few EeV [77]; sensitivities 
for 3 yrs of data taking of ANTARES [79] and of IceCube [64] (also shown for an 
E~~^ flux). Also some values of the mrf (the model rejection factor is defined as 
the average upper limit at 90% c.l. for many possible experimental outcomes over 
the expected signal) for IceCube are indicated. Notice that even if the Waxman 
and Bahcall prediction for GRBs [100] appears to be below IceCube sensitivity, the 
search for v signals in coincidence with accompanying GRBs is background free, 
resulting in a mrf < 1. On the right: 90% c.l. upper limits on diffuse v flux 
for cascades of all flavors. The AMANDA-II (197 d, indicated by 'this work') [80] 
and Baikal (268 d) [62] results are presented compared to atmospheric v fluxes, an 
AGN model (SDSS) and an upper limit extrapolated from CR UHE measurements 
(MPR). For comparison the AMANDA B-10 limits are shown multiplied by a 
factor of 3 (the underlying hypothesis is that at source the proportion between 
flavors \s, Vf, : V : Vt = \ : 1 : while at Earth, after propagation and oscillations, 
it is 1 : 1 : 1). 
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Fig. 6. On the left: fraction of e'^ measured by HEAT (circles), and AMS expected 
performances (squares) [91]. The solid line is the sum of the expected secondary 
spectrum from standard propagation models (dashed line) and the signal expected 
from annihilation of neutralinos with mass of 336 GeV. The signal has been mul- 
tiplied by a factor of 11.7 to fit the HEAT data [90]. On the right: u induced // 
fluxes induced by neutralino annihilation in the core of the Earth. Susy models are 
calculated using DarkSUSY [96] and the different symbols used indicate: models 
excluded by Eidelweiss direct search experiment [101] (green circles), models that 
would be excluded with a factor of 10 increased sensitivity (blue crosses) and those 
which would require a larger sensitivity. Experimental limits are scaled to an energy 
threshold of 1 GeV. The limits presented at the conference are indicated with the 
name of the experiment: AMANDA B-10 [95], SK [75] and Lake Baikal [62]. The 
dashed red line is the MACRO limit [97] . 



